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Abstract 

Assuming that supersymmetry is realized with parameters in the hyperbolic branch/focus point 
(HB/FP) region of the minimal supergravity (mSUGRA) model, we show that by searching for 
multijet + £vp lss events with tagged b jets the reach of experiments at the LHC may be extended 
by as much as 20% from current projections. The reason for this is that gluino decays to third gen- 
eration quarks are enhanced because the lightest neutralino has substantial higgsino components. 
Although we were motivated to perform this analysis because the HB/FP region is compatible 
with the recent determination of the relic density of cold dark matter, our considerations may well 
have a wider applicability since decays of gluinos to third generation quarks are favoured in a wide 
variety of models. 
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The WMAP collaboration has determined the cosmological density of cold dark matter 
(CDM) to be, 

VcDMh 2 = 0.11261K • (1) 

This measured value is of the same order of magnitude as the density expected from the 
production in the Big Bang of a stable, weakly interacting particle with a mass 0(100 GeV), 
assuming only that it was in thermal equilibrium at some point in the past |2|. Hence, 
the precise WMAP measurement provides a stringent constraint on all models that include 
heavy, stable weakly interacting particles. However, since the dark matter may well be made 
of several components, strictly speaking (0) only implies an upper bound on the density of 
such particles. In particular this bound applies to the stable lightest supersymmetric particle 
(LSP), frequently the lightest neutralino Z\, of .R-parity conserving supersymmetric models 
that have been the focus of much attention during the last twenty- five years [3j], and leads 
us to conclude that, 

VL^h 2 < 0.129 (2a) . (2) 

Compatibility with (j2J is possible only if the neutralinos can annihilate efficiently which, in 
turn, is possible only if one of the following holds: 

i. The lightest neutralino is hypercharge gaugino-like and annihilates via t-channel ex- 
changes of relatively light (~ 300 GeV) sfermions 

ii. The neutralino mass ~ \^la,h so that it annihilates resonantly via the exchange 
of the neutral Higgs bosons A or H in the s-channel j^J. Since these heavier Higgs 
bosons are typically quite broad, and because the neutralinos have thermal motions, 
resonant Higgs annihilation occurs over a rather large range of parameters. There is 
also a small range of parameters where resonant annihilation via the lightest scalar 



Higgs boson h leads to efficient neutralino annihilation 



iii. The neutralino mass is close to that of a charged or coloured particle; since this latter 
particle can annihilate efficiently, the neutralino density is correspondingly reduced 
as long as its interactions maintain it in thermal equilibrium with the co-annihilating 
charged particle 

iv. The parameter is small compared to the gaugino masses so that the lightest neu- 
tralino has significant higgsino components and annihilates effectively via couplings to 



electroweak gauge bosons 



v. The neutralino has significant SU(2) gaugino components, and so annihilates to 
W + W~ via large isotriplet SU{2) couplings. In this case, the lighter chargino tends 
to be close in mass to Z\, so co- annihilations may also be important jjj. 

Most SUSY analyses of the implications of the WMAP measurement have been performed 
within the framework of the mSUGRA model which, assuming radiative electroweak 
symmetry breaking, is specified by the well known parameter set, 

m , mi/2, A , tan/3, sign(/x). 

Within this framework, which has also been the paradigm for many phenomenological anal- 
yses of supersymmetry (SUSY), neutralino annihilation as in item|IJ occurs in the so-called 
bulk region with small values of m and m^; as in|nj in the A or H funnels which occur only 
if tan j3 is large; as in |m| only close to the boundary of parameter space where f\ becomes 



the LSP 



111 ], or for special values of A , where t\ becomes the LSP 12( ; and as in|Iv| in the 



hyperbolic branch/focus point (HB/FP) region with large values of tuq and modest to large 



values of /;/ J/2 |l^|. The last option v. is uol realized within mSCGR A or For that matter 
in any model with unification of gaugino masses, or in any SUSY Grand Unified Theory 
(GUT) unless the field that breaks SUSY also breaks the GUT gauge symmetry In 
the recently studied non- universal Higgs masses (NUHM) extensions of mSUGRA 15| . the 
Higgs funnel occurs for any value of tan/? while the low ]//] region is accessible even for 
relatively low values of tuq [161 ] . 

Within the mSUGRA framework that we adopt for this study, it has been shown that, 
with an integrated luminosity of 100 fb _1 , experiments at the Large Hadron Collider (LHC) 
will probe all of the bulk region, most of the Higgs funnel and, except for the largest values 

n 

of tan/?, all of the stau co-annihilation region allowed by the WMAP data [17]. The reach 
in the low \fi\ HB/FP region, however, cuts off around rrig ~ 1.6 — 1.8 TeV, where the signal 
from gluino pair production becomes rate limited. Although W% and Z 2 are relatively light 
and will be abundantly produced at the LHC, if jaI is indeed small the efficiency for the well 
studied trilepton signal from W\Z 2 production 18| is reduced, especially if we are deep in 
the HB/FP region where the mass gap between W% (or Z%) and Z\ becomes rather small and 
the daughter leptons as well as E™ 1SS from charginos and neutralinos are soft. In this case, it 
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has been shown |19j that by implementing specially designed cuts to separate the chargino 
pair production soft decay products from Standard Model (SM) background, experiments 
at an e + e~ linear collider operating at yjs = 500 — 1000 GeV will be able to probe portions 
of the HB/FP region not accessible at the LHC. Since the LHC is scheduled to commence 
operations in 2007, while a linear collider is even very optimistically at least a decade away, 
it is clearly worthwhile to explore all strategies that can potentially expand the LHC reach, 
especially in this low \fi\ region favoured by the WMAP measurements. 

An obvious option would be to re-examine the trilepton signal to see whether it is possible 
to separate it from SM background processes. 1 In this paper, however, we follow a completely 
different strategy and focus on the signal from gluino pair production (since squarks are very 
heavy). Our starting point is the observation that since the lighter higgsino-like neutralinos 
and charginos couple much more strongly to the third generation than to the first two 
generations, decays of gluino into third generation fermions will be strongly enhanced so 
that the signal may be expected to be rich in high Et 6-jets [21]. In contrast, the dominant 
SM backgrounds to multijet + E™ 1SS channels which give the largest SUSY reach at the LHC 
come from ti production, from V + jet production (V = W, Z) and from QCD processes. 
Since the latter two backgrounds are not expected to be especially rich in hard bottom quark 
jets, and because experiments at the LHC are expected to have good b tagging capability, 
we explore whether requiring the presence of tagged 6-jets in the signal allows us to probe 
portions of the hyperbolic branch that are inaccessible using the by now standard analyses 
22.0. [24 of the various multijet + E™ 1SS channels at the LHC. While 6-tagging has been 
suggested before to explore the nature of the underlying model j^j], to our knowledge it has 
never been proposed as a tool for SUSY discovery. We remark that although the 6-tagged 
signal will be rate limited, unlike the trilepton signal from W\Z 2 production, this signal will 
be relatively insensitive to how deep we are in the HB/FP region. 

We use the program ISAJET 7.69 26] with a toy calorimeter described in Ref. (3] for our 
analysis. Jets are found using a cone algorithm with a cone size AR = y/Arj 2 + A<p' 2 = 0.7. 
Clusters with Et > 40 GeV and |^(jet)| < 3 are defined to be jets. Muons (electrons) are 
classified as isolated if they have Et > 10 GeV (20 GeV) and visible activity in a cone with 



1 While we were preparing this manuscript, we learnt that the reach via the trilepton channel has recently 
been re-examined in Ref. |20j , and found to be smaller or comparable to the reach in other leptonic and 
the E™ 1SS channels even when < Mi j2 . 
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AR = 0.3 about the lepton direction smaller than Et < 5 GeV. We identify a hadronic 
cluster with Et > 40 GeV and \r)(j)\ < 1.5 as a b jet if it also has a B hadron, with 
Pt{B) > 15 GeV and \t](B)\ < 3, within a cone with AR = 0.5 of the jet axis. We take the 
tagging efficiency e& = 0.5, and assume that gluon and other quark jets can be rejected as b 
jets by a factor R b = 150 (50) if E T < 100 GeV (E T > 250 GeV) and a linear interpolation 
in between |27 1 . While we make no representation about the tagging efficiency and rejection 
against light quark and gluon jets that will be finally achieved, especially at high LHC 
luminosity and in the jetty environment of SUSY events, we felt that an exploratory study 
of just how much 6-jet tagging helps with the LHC reach would be worthwhile. 

Rather than perform time consuming scans over the WMAP favoured HP/FB regions of 
the rriQ — mi/2 planes of the mSUGRA model for several values of tan/3, we have chosen 
three diverse model lines for our analysis. For each of these model lines, we take /i > (the 
sign favoured by the result of experiment E821 at Brookhaven 
results are largely insensitive to this choice — and 



28| ) and fix A = — our 



• mi/2 = 0.295m - 477.5 GeV with tan/3 = 30 for Model Line 1, 

• mi/2 = 0.295m - 401.25 GeV with tan/? = 30 for Model Line 2, and 

• mi/2 — |o m o — 390 GeV with tan/3 = 52 for Model Line 3. 

For values of mo ^ 1500 GeV, these model lines all lie in the WMAP allowed HB/FP region 
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. The first two model lines have an 



of the mSUGRA parameter space delineated in Ref. 
intermediate value of tan /3 with Model Line 1 being deep in the HB /FP region while Model 
Line 2 closer to the periphery of the corresponding WMAP region. We choose Model Line 3 
again deep in the HB/FP region, but with a very large value of tan (3 to examine any effects 
from a very large bottom quark Yukawa coupling. We take m t = 175 GeV throughout this 
analysis. 

The branching fractions for the decays of the gluino with a mass ~ 1650 GeV (close to the 
limit that can be probed at the LHC via the usual multijet + multilepton +.E™ 1SS analyses) 
into third generation fermions are shown in Table |l] for the three model lines introduced 
above. The following features of Table |l] are worth noting. 

1. In all three cases, almost 90% of the gluino decays are to the third generation, so that 
we expect very hard top and bottom quark jets in SUSY events. 



TABLE I: The branching fractions for the decays of the gluino with m g ~ 1650 GeV to third 
generation quarks in the three WMAP allowed model lines introduced in the text, together with 
the corresponding values of the hypercharge gaugino mass parameter Mi, and the superpotential 
Higgs mass parameter [i. 

Model Line 1 Model Line 2 Model Line 3 



m~ g (GeV) 


1637 


1665 


1665 


fi (GeV) 


181 


335 


164 


Mi (GeV) 


269 


278 


275 


g -> W{tb + Wfib 


37% 


34% 


39% 


g -> W 2 ~tb + W 2 + ib 


8% 


11% 


8% 


g Zitt + 


15% 


5% 


17% 


5 -» Z 2 tt + Z 2 bb 


19% 


16% 


20% 


g -» Z 3 tt + Z 3 bb 


6% 


17% 


4% 


g —* z^tt + 


4% 


5% 


4% 



2. In Model Lines 1 and 3 that are deep in the HB/FP region, the gluino mainly decays 
to the higgsino-like lighter chargino and the two lightest neutralinos; in Model Line 2, 
Z\ is dominantly the hypercharge gaugino, so that gluino decays to Z 2 and to Z 3 are 
favoured. 

3. The main difference due to the large tan/? value for Model Line 3 is the increased 
branching ratio for the decays g — > bbZi relative to g — > ttZi. Although we have not 
separated these out in the Table, we have checked that while the direct decays to 
bottom comprise just about 10% of all gluino decays to neutralinos for Model Lines 
1 and 2, these decays constitute about a third of all gluino to neutralino decays for 
Model Line 3. This is, of course, due to the increased Yukawa coupling of the bottom 



quark (29 ]. 



4. In each of these cases, we see that decays of the gluino to the wino-like charginos 
and neutralinos have relatively small branching fractions despite their large SU(2) 
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gauge couplings. This is because decays mediated by lighter third generation squarks 
that have large Yukawa couplings dominate because these are dynamically as well as 
kinematically favoured. 

The main message of this Table is that in the WMAP favoured HB/FP regions of the 
mSUGRA model, decays to third generation quarks dominate gluino decays. Moreover, 
although detailed decay patterns depend on both tan/? and the value of ///Mi ( i.e. on 
how deep we are in the HB/FP region), the total branching fraction for decays to third 
generation quarks is relatively insensitive to these details. Motivated by these observations 
we begin our examination of the inclusive 6-jet signal for each of the model lines introduced 
above. 

The major SM backgrounds to the multijet plus E^ 1SS signal, with or without 6-jets, 
come from W or Z + jet production, from ti production and from QCD production of 
light quarks and gluons. In the last case, the i?™ lss arises from neutrinos from c and b 
quarks, from showering of W, Z bosons and their subsequent decays to neutrinos, and from 
energy mismeasurement. The backgrounds from these sources are shown in Table ITT1 for two 
representative choices of cuts discussed below. Here St is the transverse sphericity, m e g 
is the scalar sum of the Et of the four hardest jets and E^ 1 ^, A0 is the azimuthal angle 
between i?™ lss and the hardest jet, and is the azimuthal opening angle between the two 
hardest b jets in events with Nj, > 2. Since we do not require multileptons in our analysis, 
backgrounds from WW, WZ and ZZ as well as three vector bosons processes are expected 
to be much smaller J^]. 

We see that for both sets of cuts the SM background is dominated by QCD. This differs 
from earlier results where it is argued that the QCD background can be reduced to negligible 
levels by requiring the signal to be sufficiently hard. Indeed, to track the reason for this 
difference, we have done a very high statistics analysis of the QCD background that was not 



possible in Ref. 22]. Specifically, we took particular care to divide the event generation into a 
large number (50 bins for QCD, fewer for other backgrounds) of finely spaced hard scattering 
Pt bins, especially for lower values of the hard scattering pt where the event weights are 
very large, even if the efficiency for passing the cuts is low. 2 We have checked that our 



2 If, for any set of cuts, we find no events in our simulation of a particular background, we set this background 
cross section to a value corresponding to the single event level in the smallest weight bin in our simulation. 
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TABLE II: Cross sections for the dominant SM backgrounds to SUSY processes at the CERN LHC 
for two representative choices of cuts, together with signal cross sections for three cases along the 
model lines introduced in the text. 



Variables 




Cut 1 (cut 2) 








4(5) 




Ef iss (GeV) > 




400 (500) 




Eft (GeV) > 




400 (400) 




Eft (GeV) > 




250 (250) 




Eft (GeV) > 




150 (175) 




Eft (GeV) > 




100 (125) 




m cS (GeV) > 




1500 (2250) 




S T > 




0.0 (0.0) 




Acj) < 




180° (140°) 








180° (180°) 




N b > 





1 


2 


Source 




<r (fb) 




QCD 


27.71 (0.636) 


8.95 (0.112) 


1.72 (0.009) 


ti 


5.17 (0.043) 


3.42 (0.034) 


1.10 (0.005) 


W + jets 


13.60 (0.462) 


1.42 (0.045) 


0.13 (0.002) 


Z + jets 


5.68 (0.180) 


0.65 (0.018) 


0.05 (0.001) 


Total 


52.16 (1.32) 


14.44 (0.209) 


3.00 (0.017) 


Case 1: m- g = 1054 GeV 


12.1 (0.667) 


9.25 (0.537) 


4.58 (0.290) 


Case 2: m- g = 1436 GeV 


3.81 (0.541) 


3.10 (0.47) 


1.68 (0.24) 


Case 3: m- g = 1705 GeV 


1.11 (0.230) 


0.90 (0.185) 


0.50 (0.106) 



backgrounds levels from W + jets, Z + jets and ti processes are in agreement with those in 
Ref. [22j| and attribute the difference in the QCD background to statistics of the background 
simulation. 3 Since QCD typically contributes about half the background in Table |HJ we 

3 In passing, we mention that we also found a significant QCD contribution to the background in the 



s 



may expect a small degradation of the LHC reach from these earlier projections. 4 

Also shown in Table |H] is the SUSY signal for three points along model lines. In our 
analysis, we consider a signal to be observable with a given integrated luminosity if, (a) its 
statistical significance N s /y/T^B > 5, (b) N S /N B > 0.25, and (c) N s > 10. We see that 
Case 1 with rrig ~ 1 TeV which by early analyses should be easily observable at the LHC 
is also observable using cuts 1 with an integrated luminosity of just 10 fb _1 , and moreover, 
the significance of the signal improves with increasing b multiplicity. For Case 2, the SM 
background with the softer cuts 1 is too large except in the 2b channel, where the signal 
is observable for an integrated luminosity > 26 fb . However with the harder cuts of set 
2 the signal, though unobservable without 6-tagging, should be observable in the lb (2b) 
channel for an integrated luminosity exceeding 26 fb -1 (41 fb -1 ). For Case 3, both 6-tagging 
capability and an integrated luminosity of at least 72 fb -1 are essential for the observability 
of the signal. It is clear that 6-tagging improves the reach of the LHC for points in the 
HB/FP region. 

To quantify the improvement 6-tagging makes to the capabilities of the LHC for the 
detection of SUSY, we have re-computed the reach of the LHC for each of the three model 
lines introduced above. Towards this end, for every mSUGRA parameter point that we 
examined, we generated a set of SUSY events using ISAJET. We also generated large samples 
of SM background events. We passed these events through the toy calorimeter simulation 
mentioned previously, and then analysed both the signal and the background for the entire 
set of cuts (5 x 5 x 6 x 5 x 3 4 = 60750 choices in all) listed in Table IIII1 We regard the 
signal as observable if it satisfies the observability criteria (a)-(c) listed above for at least 

multijets + \l + E™ lss channel. This is in contrast to the result in Ref. [30( which was obtained using 
PYTHIA. We attribute this difference to the fact that PYTHIA does not include showering of W and 
Z bosons in QCD events. While this may result in some double counting of "Drell-Yan" W and Z 
production, we note that showering of vector bosons from the final state quarks will, in general tend to 

populate a different region of phase space. 
4 A more significant background issue may be that the showering algorithms appear to obtain a significantly 

softer distribution of the variable rru; (introduced below) relative to evaluations using exact multijet plus 
W, Z production matrix elements [31|. We have nothing to say about this, except that requiring additional 
6's (and in the case of W, also requiring a transverse mass cut) will reduce this background considerably. 
At the very least, our analysis will indicate the extent to which the presence of tagged b's increases the 
LHC reach. We may add that the E^ analysis of Ref.|22| that requires E^ lss together with just two 
(rather than four) additional hard jets, with just one of these jets from showering, may be more robust 
to these matrix element corrections. 
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TABLE III: The set of cuts used to optimize the SUSY signal at the LHC along the hyperbolic 
branch of the mSUGRA model. 



Variable 


Values 


Nj> 


4, 5, 6, 7, 8 


Ef iss (GeV) > 


400, 500, 600, 700, 800 


(E T \ E T 2 , E J T 3 , E J T 4 ) (GeV) > 


(400, 250, 150, 100), 




(400, 250, 175, 125), 




(400, 250, 200, 150), 




(500, 350, 150, 100), 




(500, 350, 175, 125), 




(500, 350, 200, 150) 


m cS (GeV) > 


[E™ ss + Eii + 200 x n, n = 0, 1, 2, 3, 4 


S T > 


0.0, 0.1, 0.2 


Acj) < 


180°, 160°, 140° 


A<p b < 


180°, 160°, 140° 


N b > 


0, 1, 2 



one choice of cuts. 

Notice that the cuts in this Table are much harder than those used in a recent analysis of 
the LHC SUSY signal [l^. This is because, unlike in Ref. jl^ where the cuts were designed 
to extract the signal for a wide range of squark and gluino masses, here we focus on the 
optimization of the signal in the portion of the HB/FP region with heavy gluinos where 
the previous strategy fails. Note also that cut 1 in Table [TT1 corresponds to the softest of 
these cuts. Since, as we saw earlier, the signal for Case 1 (with rru ~ 1 TeV) is comfortably 
observable with the present set of cuts as well as those in Ref. [17[, there is no danger that 
there will be a gap in the HB/FP region of parameter space where the signal is unobservable 
with either strategy. 

Our results for the LHC reach are shown in Fig. ^ where we plot the largest statistical 
significance of the signal as we run over the cuts in Table ITTT1 for a) Model Line 1, b) Model 
Line 2, and c) Model Line 3. The solid curves, from lowest to highest, denote this maximum 
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FIG. 1: The maximum value of the statistical significance for the multijet + Ej^ ms SUSY signal at 
the LHC as we run over the cuts listed in Table ITTTI for the three model lines introduced in the text, 
requiring in addition that N$ > 10 events and that Ng > 0.25Nb- The solid lines from bottom 
to top denote this statistical significance without any b tagging requirement, with at least one 
tagged b-jet, and with at least two tagged 6-jets respectively, assuming an integrated luminosity of 
100 fb _1 and a tagging efficiency of 50%. The dotted lines show the corresponding result for an 
integrated luminosity of 300 fb _1 . 

statistical significance without any b tagging requirement, requiring > 1 tagged 6-jet, and 
> 2 tagged 6-jets respectively, for an integrated luminosity of 100 ftr 1 , while the dotted 
curves show the corresponding results for 300 fb _1 of integrated luminosity that may be 
expected from three years of LHC operation with the high design luminosity. While, without 



6-tagging, our reach for gluinos is ~ 200 GeV smaller than earlier projections 17, |23j,|2j, 
presumably because of differences in the background levels discussed above, we see that 
6-tagging will improve the mass reach of gluinos by 15-20%, provided that LHC experiments 
can accumulate an integrated luminosity of 100-300 fb _1 and that 6-tagging with an efficiency 
of ~ 50% remains possible even in the high luminosity environment. 
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A few remarks appear to be in order at this point: 

• The statistical significance in Fig. ^ is not significantly improved if the 6-tagging 
efficiency improves to 60%. The reason is that before tagging the signal typically 
contains (on average) 3-4 b quark jets while the background typically contains (at 
most) just two b quark jets. As a result, the increased efficiency enhances the b- 
tagged background more than the signal, and the statistical significance is essentially 
unchanged. 

• We have also checked that with a 6-tagging efficiency of 50% and an integrated lumi- 
nosity of 100 fb _1 , the signal with > 3 tagged 6-jets is rate limited and no increase in 
the reach is obtained from that shown in the Figure. 5 

• The search strategy proposed here does not use lepton information at all. We have 
checked that a transverse mass cut Mt(£, E™ 1SS ) > 100 GeV on events with at least 
one isolated lepton does not increase the significance of the signal because the fraction 
of signal events (after our cuts) with an isolated lepton is not especially large. 

• Although we have not shown this explicitly, we have checked that requiring the pres- 
ence of additional isolated leptons does not lead to an increase in the reach relative to 
the > 26 channel. 

In summary, we have shown that in the HB/FP region of the mSUGRA model the reach 
of the LHC as measured in terms of gluino masses may be increased by 15-20% by requiring 
the presence of hard, tagged 6-jets in SUSY events. While we were mainly motivated in our 
investigation by the fact that this part of parameter space is one of the regions compatible 
with the WMAP data, our considerations may have wider applicability since decays of 
heavy gluinos to third generation fermions are favoured in all models with common masses 
for sfermions with the same gauge quantum numbers. This is in part because the large top 
Yukawa coupling and, if tan (3 is large, also the bottom quark Yukawa coupling, cause the 
third generation squarks to be lighter than their siblings in the first two generations, and 

5 It is possible that a slightly increased reach may be obtained if the tagging efficiency is significantly larger 
than 50% or if the integrated luminosity is considerably higher than 100 fb _1 . It would, however, be 
necessary to evaluate backgrounds from 46, 4i and ttbb production processes before a definitive conclusion 
can be made. 
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in part because of new contributions to gluino decay amplitudes from these large Yukawa 
couplings (s^l. Specifically, we may expect that signals with tagged 6-jets may also be 
useful in models with an inverted squark mass hierarchy 32[, in models with unification 
of Yukawa couplings (because they require large tan/3), and possibly also in models with 
non-universal Higgs mass parameters that have recently been re-examined in light of the 
WMAP data Q]. Since the CMS and ATLAS experiments are expected to ultimately have 
good 6-tagging capability, we urge that it be utilised to maximize the SUSY reach of the 
LHC. 
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